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Abstract 

We investigate the nonlinear response of an anharmonic monostable quantum me- 
chanical resonator to strong external periodic driving. The driving thereby induces 
an effective bistability in which resonant tunneling can be identified. Within the 
framework of a Floquet analysis, an effective Floquet-Born-Markovian master equa- 
tion with time-independent coefficients can be established which can be solved 
straightforwardly. Various effects including resonant tunneling and multi-photon 
transitions will be described. Our model finds applications in nano-electromechanical 
devices such as vibrating suspended nano-wires as well as in non-destructive read- 
out procedures for superconducting quantum bits involving the nonlinear response 
of the read-out SQUID. 
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1 Introduction 



Classical nonlinear systems subjected to strong periodic external driving of- 
ten have several stable stationary states for which the amplitudes and phases 
of the forced vibrations differ in size [1,2,3]. One of the simplest theoretical 
models which show the coexistence of two stable states induced by external 
driving is the well-known classical Duffing oscillator. An anharmonic statically 
monostable potential can be driven into a dynamically bistable regime showing 
various interesting features of non-linear response [1,2,3], such as hysteresis, 
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period doubling, and thermal activation when finite temperatures are consid- 
ered. The external driving field with frequency uo induces an effective dynamic 
bistability which is manifest by the non-monotonous dependence of the am- 
plitude A of the stationary vibrations for varying id. For the classical system 
where all potential energies are allowed, this response curve A(u>) is smooth 
showing only two points of bifurcation for the related bistability. If the control 
parameter u is additionally varied adiabatically, hysteretical jumps between 
the two stable states occur. If additional thermal noise is added to the system, 
the regime of bistability shrinks due to thermal escape of the metastable state. 

The main subject of this work is to investigate the corresponding driven quan- 
tum mechanical system. The presence of time-dependent driving typically adds 
several interesting features to the properties of the time-independent quantum 
system, see for instance Refs. [4,5]. In this work, the focus is laid on the non- 
linear response of the driven quantum mechanical anharmonic oscillator in the 
presence of an Ohmic heat bath. We show that the nonlinear response curve 
A(uj) displays beyond its characteristic shape additional quantum mechani- 
cal resonances which are related to multi-photon absorptions. Additionally, 
we show that there exists a separation of time-scales indicating that the sta- 
tionary state is reached by quantum tunneling from the dynamically induced 
metastable to the globally stable state. In fact, by tuning the control param- 
eter u, several resonant tunneling transitions can be identified as resonances 
in the corresponding tunneling rate. 

The paper is organized in the following way: In Section 2, the system Hamilto- 
nian is introduced. The latter is periodic in time, which allows the application 
of Floquet theory. Since we are interested in the stationary state in presence 
of (weak) dissipation, we introduce also a set of harmonic oscillators repre- 
senting an Ohmic heat bath. An efficient way to determine the dynamics of 
the system is the use of a Born-Markovian master equation in the Floquet pic- 
ture (Section 3) yielding a simple master equation with time-independent rate 
coefficients. After straight forward diagonalization, the stationary oscillation 
amplitude A and the phase tp follow. In Sections 4 and 5, the nonlinear re- 
sponse of A and ip depending on the various parameters is studied in detail. In 
Section 6, resonant tunneling is investigated. Finally, in Section 7, we discuss 
the applicability of our model to experimental systems before we conclude. 



2 The quantum Duffing oscillator 

The Hamiltonian for the driven anharmonic oscillator has the form 

H s (t) = 7^ + -^x 2 + -x 4 + xf cos(cut) . (1) 
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Here, m and ujq are the mass and the harmonic frequency of the resonator, 
respectively, while a gives the strength of the nonlinearity. We focus on the 
case a > of hard non-linearities, where the undriven potential is monostable. 
The external driving is characterized by the amplitude / and the frequency 
uj. As it will become clear below, the driving induces an effective bistability 
in which quantum tunneling can be identified. 

We include the effect of the environment by a bath of harmonic oscillators 
which are bilinearly coupled to the system with the coupling constants Cj [6]. 
The Hamiltonian for the bath and the coupling to the system is given by its 
standard form 

We focus on the generic case of an Ohmic bath with the spectral density 

7T C 2 

J(uj) = -J2 -Uj)= m 7 we- u/Wc , (3) 

j m i UJ i 

with damping constant 7 and cut-off frequency uo c . The total Hamiltonian is 
H(t) = H s (t) + H B . 

To proceed, we scale H(t) to dimensionless quantities such that the energies 
are in units of Hloq while the lengths are scaled in units of x = \J^j^- Put 
differently, we formally set m — H — uq — 1. The nonlinearity parameter a is 
scaled in units of a = hw /xQ, while the driving amplitude is given in units of 
/ = hu /x . Moreover, we scale temperature in units of T = fiw /k B while 
the damping strengths are measured with respect to c^o- 

To investigate the dynamical behavior of the driven resonator, it is convenient 
to use the periodicity of H s (t) with respect to time and switch to the Floquet 
picture [7], the later being equivalent to a transformation to the rotating frame. 
The Floquet or quasi-energies e a follow from the solution of the eigenvalue 
equation 

H s {t)-ih^ \<f> a (t)) = e a \<j) a (t)} , (4) 

with the Floquet states \(f> a (t)). The quasienergies e a are defined up to a mul- 
tiple integer of fko, i.e., the state \<f>£\t)) = e mu)t \(j) a (t)) is also an eigenstate 
of the Floquet Hamiltonian, but with the eigenvalue £ Q>n = e a + nhio. This 
feature prevents us from a global ordering of the quasienergies which, how- 
ever, can be achieved with the mean energies obtained after averaging over 
one driving period, i.e., 

), (5) 
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with the Fourier components of the Floquet states [7] 

\c a , n ) = — dte tnut \Mt))- (6) 

Z7T JO 

3 Dynamics of the quantum Duffing oscillator 

3.1 Floquet-Born-Markovian master equation 

The dynamics of the resonator in the regime of weak coupling to the bath can 
be efficiently described by a Born-Markovian master equation in the Floquet 
picture [7] for the elements p a p{t) = (<f> a (t)\p\ fipit)) °f t ne reduced density op- 
erator p after the harmonic bath has been traced out in the usual way. For weak 
damping, the dissipative influence of the bath is relevant only on a time scale 
much larger than the driving period T w = 2-k/uo. Thus, the time-dependent 
coefficients which are periodic in time with period T w can safely be replaced 
by their average over one driving period (moderate rotating wave approxi- 
mation [7]). This yields a simplified master equation with time-independent 
coefficients which reads 

Pap(t) ^Y, M ocp,a'p'Pa'p'(t) (7) 
a'P' 

with 

■MaP,a'P' = — ^ [£ a — e p) baa'Spp' + £ a p,a'P> ■ (8) 

The first term on the r.h.s. describes the coherent time evolution of the pure 
system while the second term contains the transition rates describing the 
influence of the dissipative bath. It reads [7] 

CaP,a'P> = X! (Naa',n + ^PP',n) X aa i , n X pi p _ n — Spp* ^ ^ P" cx' ,n^- a,p" -n^- P" o! ,n 
n p", n 

— $aa' ^2 N a " P' ,nXp> a " ,-n^a" P,n ■ (9) 
a" ,n 

with the coefficients 



TWyE 1 

N aP , n = N(e a -e p + nfuv), N(e) = - 2 ^ /kBT _ - , 



X a p, n = — / dte- mu,t (cf) a (t)\x\Mt))- (10) 

Z7T JO 

The operator Ai in Eq. (8) is a M 2 x M 2 — matrix, when the total Hilbert 
space of the anharmonic resonator has been truncated to the M— dimensional 



4 



subspace. For practical purposes, we set M — 12 throughout this work. Note 
that we have confirmed convergence with respect to M for all results shown 
below. With that, M. can be readily diagonalized numerically by standard 
means. This can be formalized in terms of the diagonalization transformation 
S by the eigenvalue equation 

(S 1 )aP, f MuM^ fl > u >S^ u > :a /[3i = A a p5 ata >5[3 t i3> . (11) 

/ii'j/iV 

Here, A a/ 3 denote the eigenvalues of the operator M.. Since the master equation 
(7) conserves the trace of p, there is an eigenvalue = which character- 
izes the stationary solution p^. The remaining eigenvalues all have a negative 
real part leading to a decay of the corresponding mode with time. Due to the 
structure of the master equation (7), there exist two classes of eigenvalues: (i) 
the first class having an imaginary part of zero consists of individual eigen- 
values (associated to relaxation), and (ii) the second class having non-zero 
imaginary parts (associated to dephasing) consists of pairs of complex conju- 
gated eigenvalues. The eigenvalues can be ordered according to the size of the 
absolute value of their real parts. As it turns out, there is one pair of complex 
conjugated eigenvalues, i.e., A^g in = A^* min with the smallest non-zero ab- 
solute value of the real part. Moreover, we find that that pair of eigenvalues is 
clearly separated in size from the remaining ones with respect to the real parts 
indicating a separation of time scales. That pair of eigenvalues is responsible 
for the tunneling dynamics at long times, as will be discussed in Section 6. 
Finally, the solution of the master equation (7) can be formally written as 

P^(t) = Pap+ E S a ^(S- 1 ) llu ^,e A ^ t p fl , u ,(t = 0). (12) 
For convenience, we have chosen the initial time t — 0. 

3.2 Observable: Stationary oscillation amplitude 

We are interested in calculating the asymptotic expectation value (x) ttSlS of 
the position operator. This is the quantity which can directly be compared to 
its classical counterpart being the solution of the classical Duffing equation. 
It reads 

(x) t , as = J2o n e muJt (13) 

n 

with the coefficients 

C n = A n e i<Pn = Paf3 X f3a,n ■ (14) 
a/3 

Since the oscillator is driven by a cosine-shaped external force, it will oscil- 
late in the asymptotic limit also with a phase-shifted cosine if the excitation 
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frequency u remains close to uq. Higher harmonics could be generated in this 
nonlinear system but they are not in the focus of our interest in this work. 
Thus, only the terms n — ±1 contribute in the Fourier expansion and we 
obtain 

( x )t,as — A COs(wt + <f) , OJ ~ UJq (15) 

with the amplitude A = 2A 1 and the phase ip = ipx of the first harmonic of 
the Fourier expansion. These two quantities are used to study the nonlinear 
response of the anharmonic resonator in the stationary long-time limit. The 
short time dynamics of such a type of master equation is an interesting issue 
by itself since it is related to the question of complete positivity and to the 
question of factorizing initial conditions [8,9]. Moreover, we note that the 
master equation (7) is valid only in the case of weak system-bath coupling. For 
the opposite limit of strong coupling (quantum Smoluchowski limit), different 
techniques [10,11,12,13,14] have to be applied. 

3.3 Classical Duffing oscillator 

The corresponding classical oscillator (at zero temperature) is the well-known 
Duffing oscillator [1,2,3]. It shows a rich variety of features including regular 
and chaotic motion. In this work, we focus on the parameter regime where 
only regular motion occurs. The nonlinear response of its amplitude A can be 
calculated perturbatively [1]. One obtains the response A(uj) as the solution 
of the equation 

3 a P 7 2 V /2 ritt 

Its characteristic form is shown in the inset of Fig. 1. For weak driving 
strengths, the response as a function of the driving frequency u has the well- 
known form of the harmonic oscillator with the maximum at uj = ujq. For 
increasing driving strength, the resonance grows and bends away from the 
uj = cj -axis towards larger frequencies (since a > 0). The locus of the maxi- 
mal amplitudes is given by the parabola [1] uj — uj = ^^^A 2 , which is often 
called the backbone curve. Most importantly, a bistability develops with two 
adjacent stable branches and one intermediate unstable branch. This bistabil- 
ity is connected with a hysteretical jump phenomenon which can be probed if 
the driving frequency u is adiabatically increased or decreased. The hysteresis 
is maximal for zero temperature. At finite temperature, it is reduced since the 
particle can escape from the metastable local minimum to the adjacent global 
minimum via thermal hopping [15,16,17] before the deterministic switching 
point is reached [18,19]. Note that for stronger driving amplitudes, also bifur- 
cations and period doublings can occur [1,2,3] which we do not address in the 
present work. 
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The nonlinear response of the phase tp can be determined perturbatively in a 
similar way. One obtains [1] 



tp — — arctan 



(17) 



2 ((a; - cu )A 
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where A is the solution of Eq. (16). The curve <p(uj) also has two stable branches 
with an unstable intermediate branch and displays similar hysteretical jump 
phenomena as the amplitude response. 

As we will show in the following, the nonlinear response is qualitatively dif- 
ferent in the corresponding quantum system. The discrete quasi-energy spec- 
trum allows for multi-photon excitations which yield discrete resonances in 
the amplitude response profile. Moreover, the dynamically generated bistabil- 
ity allows for an escape of the system out of the metastable state via resonant 
quantum tunneling. This generates characteristic resonances in the tunneling 
rate when the external frequency w, which plays the role of a control param- 
eter, is tuned. 



4 Amplitude response 

A typical response profile for the amplitude A(u) is shown in Fig. 1. The 
shoulder-like shape is a remnant of the classical form of the response which 
is indicated by the dashed line. In the quantum case, clear resonances can be 
observed at particular values of the driving frequency. 

The resonances can be understood as discrete multi-photon transitions occur- 
ring when the iV-th multiple of the field quantum fkv equals the correspond- 
ing energy gap in the anharmonic oscillator. By inspection of the associated 
quasienergy spectrum shown in Fig. 2, one can see that the distinct reso- 
nances occur at multiple degenerate avoided level crossings of the quasienergy 
level En of the A-th Floquet state with the quasienergy level So of the Flo- 
quet groundstate. In physical terms, multi-photon excitations occur in the 
resonator and the corresponding populated Floquet state dominates then the 
position expectation value (x) with its large value of the amplitude, see Eqs. 



The width of the A-th resonance is determined by the minimal splitting at 
the avoided quasienergy level crossing which is equal to the A-photon Rabi 
frequency. The latter can be evaluated perturbatively for weak driving ampli- 
tudes [20,21] upon applying a rotating- wave approximation (RWA). For the 
resulting Hamiltonian, one finds that to lowest order in f /{uj — cun), the A- 
photon Rabi frequency at decreases exponentially with A implying that 



(13,14). 
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the resonances are sharper for larger N. For small frequencies, the broad peaks 
overlap strongly and lead to a shoulder-like profile which is similar to the clas- 
sical result (dashed line in Fig. 1). Although the RWA yields a qualitatively 
correct picture, a quantitative comparison shows noticeable deviations from 
the results for the full anharmonic resonator and is not pursued further in 
this work. Note that a similar system has been investigated in Ref. [22] in the 
context of a dispersive optical bistability. There, a nonlinear Hamiltonian has 
been derived from the Duffing oscillator using a RWA, similar to that used in 
Ref. [23,20]. A rather involved matrix continued- fraction method also revealed 
a bistability. However, the numerical procedure to obtain those results is rather 
cumbersome. In contrast, our approach is numerically straightforward since it 
only involves a simple numerical diagonalization of a matrix. 

4-1 Tuning system parameters 

Since the location of the multi-photon resonances is determined by the system 
Hamiltonian, it is interesting to see how they depend on varying the system 
parameters a and /. 

4-1.1 Varying the nonlinearity coefficient a 

The spectrum of the system and therefore, the Floquet spectrum, depends 
strongly on the nonlinearity coefficient a. Increasing a increases the energy 
gaps between the succeeding eigenstates. Thus, the multiphoton absorption 
processes occur at larger frequencies of the driving field. This behavior is shown 
in Fig. 3. The increase of a leads in general to a shift of the response curve 
A(uj) towards larger frequencies u. Moreover, we observe that the height of the 
iV-photon peak decreases for increasing a. This observation is qualitatively in 
agreement with a perturbative treatment for weak driving within a rotating 
wave approximation [20]. 

4-1-2 Varying the driving amplitude f : Multiphoton antiresonance 

The Floquet spectrum of the uncoupled system depends on the value of the 
driving amplitude /. A perturbative analysis [20] shows that the iV-photon 
Rabi frequency fl^N which is given by the minimal splitting of the quasienergy 
at the avoided quasienergy level crossing depends crucially on iV and /. This, 
in turn, determines the behavior of the iV-photon resonance. The detailed 
results of the dependence on the / are shown in Fig. 4. The N = 5-photon 
peak grows for growing /. Most interestingly, the N = 6-photon peak displays 
a nontrivial behavior. For weak driving, a 6-photon antiresonance, i.e., a dip, 
develops. When the field amplitude is increased the antiresonance turns into a 
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true resonance which grows further for growing /. The particular dependence 
of the local extrema are shown in the inset of Fig. 4. As follows from Eqs. 
(13) and (14), the amplitude A is determined by the sum c 1 = J2 a p P'a^Pa,!- 
Thus, the product of the weights p^J (which also contains the influence of 
the bath) and the matrix elements Xp a! i (which is a property of the coherent 
driven system alone) determines the full shape of the amplitude. Notice that 
p^p is in general not diagonal. This is different from the system considered in 
Ref. [20] where only a coherent system without bath has been investigated. 

By varying the system parameters, we have modified implicitly both p^ and 
Xp a> -]_. In the following section, only bath parameters and thus p^ will be 
modified, leaving Xp a> i unchanged. 

4-2 Dependence on bath parameters 

The crossover from antiresonant to resonant behavior occurs when an in- 
creased driving amplitude / increases the population of higher-lying Floquet 
states. In this subsection, we investigate the role of the dissipative environment 
by tuning the bath temperature and the damping constant. 

4-2.1 Varying temperature T 

Since increasing temperature T leads also to an increased population of other 
Floquet modes, the transition from antiresonant to resonant behavior can also 
be expected to occur for growing T, at least in a certain temperature regime. 
The result is shown in Fig. 5 for the 5-photon and 6-photon-resonance. The 
resonant peak for N = 5 grows and broadens if temperature is increased. In 
contrast, the antiresonant dip for N = 6 shrinks for growing T and turns 
into a resonant peak. This peculiar behavior can be interpreted as thermally 
assisted cross-over from antiresonant to resonant behavior. The inset in Fig. 
5 shows the local extrema for N = 5 and N = 6. For even larger temperature 
k B T fkv , the characteristic peak structure is completely smeared out due 
to thermal broadening, see also Fig. 2, lower panel. Then, the peaks overlap 
and the dynamic bistability is smeared out by thermal transitions between the 
two (meta-)stable states. 

4-2.2 Varying the damping constant 7 

The results for different damping constants 7 are shown in Fig. 6. The five- 
photon resonance decreases when 7 is increased from 7 = O.OOlc^o to 7 = 
O.Olc^o. As also shown in the inset of Fig. 6, the peak maximum decreases 
monotonously for N = 5. For the six-photon (anti-)resonance, we find a dif- 
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ferent behavior. For weak damping, a sharp resonance occurs which is turned 
into an antiresonance for larger damping. The non-monotonous dependence 
of the six-photon resonance is also shown in the inset of Fig. 6. As it is the 
case for very large temperature (see above), the influence of a strong coupling 
to the bath gradually smears out and finally destroys the resonances. 



5 Phase response 

Next, we address the nonlinear response of the phase (p. The corresponding 
classical Duffing oscillator shows an interesting nonlinear phase response <p(u)) 
including two stable branches [1]. The characteristic multi-photon transitions 
in the quantum version of the Duffing oscillator also shows up in the phase 
response profile. The results for different driving strengths / are shown in 
Fig. 7 for the regime where the five- and six-photon transition occur. For 
uj < u a phase shift ip = —it is found. For u ^> u , the phase shift vanishes, 
(p — 0. In the intermediate region u > ujq, the multi-photon resonances induce 
also an antiresonance in the phase shift. The five-photon antiresonance is 
enhanced and broadened if the driving strength / is increased from / = 0.09/ 
to / = O.l/o. Increasing the driving further to / = 0.1 1/ wipes out the 
antiresonance completely and the transition ip = —n to ip — is shifted to 
larger values of cu, where this development is repeated at the succeeding six- 
photon resonance. 

This behavior is also found when temperature is varied, see Fig. 8. Increasing 
T leads to a suppression of the transition of the phase shift ip = —tt —></? = 0. 
In contrast, increasing damping favors this transition at lower frequencies u>, 
see Fig. 9. 



6 Resonant tunneling in the driving induced bistability 

The dynamic bistability of the steady state of the classical Duffing oscillator 
does not survive in the quantum system. The reason is that the system will 
escape the metastable state asymptotically via tunneling, similar to the case 
of the driven double-well potential [24]. Note also that, as a consequence, the 
hysteretical behavior is suppressed if the control parameter uj is varied truly 
adiabatically. 

Nevertheless, signatures of the dynamic bistability and tunneling can be found 
if we consider how the steady state is reached. For this, we show in Fig. 
10 the time evolution of the amplitude A (local maxima of the vibrations) 
starting with the ground state of the undriven oscillator as the initial state. We 
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observe fast oscillations at short times. They decay on a time scale 7 _1 which 
reflects "intrawell" relaxation in the metastable state. Then, starting from a 
metastable state at intermediate times, a slow exponential decay towards the 
asymptotically globally stable state can be observed. This separation of time 
scales is a clear indication of tunneling from the meta- to the globally stable 
in the dynamic bistability. 

The decay rate T for this slow process (tunneling rate) is determined by the 

(1 2) 

absolute value of the real part of the eigenvalue Aig 4in °f the operator M. in 
Eq. (8), i.e., V = | Re A^ 2 ^ in | . Results for the tunneling rate as function of the 
control parameter uj are shown in Fig. 1 1 for two different damping constants 
7. Most importantly, the tunneling rate shows resonances at the same values 
of the frequencies where the avoided crossings of the quasienergy levels occur 
(see dashed vertical lines). The peaks in Y indicate resonant tunneling [24] 
from the meta- to the globally stable state both of which are dynamically 
induced. Note the analogy to resonant tunneling in a static double-well po- 
tential [24]. The role of the eigenenergies in the static case is now played by 
the quasienergies e a determining the coherent dynamics, see Eq. (7). In both 
cases, the avoided (quasi-)energy level crossings are the origin of resonant tun- 
neling. Nevertheless, the incoherent part of Eq. (7) is crucial to observe the 
resonant tunneling in this driving induced bistability. 

Furthermore, it is interesting to note that the resonant tunneling is enhanced 
if the coupling to the bath is increased from 7 = 0.001cj to 7 = 0.005u;o 
(bath assisted resonant tunneling). We have also calculated the dependence 
of this phenomenon on temperature but we found a weak dependence in the 
interesting low-temperature regime (not shown). 



7 Discussion and conclusions 

The quantum Duffing oscillator is a generic theoretical model which finds sev- 
eral applications in experimental systems. For instance, a suspended nanome- 
chanical beam [25,26] which is excited to transverse vibrations behaves as 
a damped anharmonic resonator. Several experimental groups have observed 
the behavior described by the classical Duffing oscillator [27,28,29,30]. The 
transition to the quantum regime is currently in the focus of intense research 
[31,32,33,34,35,36]. Once, such kind of truly quantum-'mechanical' systems 
on the nanoscale have been shown to exist, macroscopic quantum effects [33] 
should be readily observable. A second class of experimental systems ad- 
dresses the resonant non-destructive read-out of a persistent current qubit 
[37,38,39,40]. In contrast to the conventional switching current measurement 
that generates unwanted quasi-particles when the dc-SQUID (acting as the 



11 



qubit detector) switches to the voltage state, this technique keeps the SQUID 
biased along the supercurrent branch during the measurement. Thereby, the 
Josephson plasma resonance of the SQUID depends on the inductive cou- 
pling of the SQUID to the qubit. Measuring the plasma resonance allows to 
non-destructively read out the qubit state. The application of this read-out 
technique in the nonlinear regime of the SQUID could allow for an improved 
sensitivity as well as its potential use as a nonlinear amplifier. Finally, we note 
that there exists a wide parameter regime (typically strong driving and/or big 
nonlinearities) where a quantum chaotic behavior of the system with many 
interesting features can occur [41,42]. A detailed study of this regime goes 
beyond the scope of this work. 

To summarize, we have investigated the nonlinear response of the amplitude 
as well as of the phase of the oscillations of a driven damped anharmonic 
resonator (quantum Duffing oscillator). The use of an efficient Floquet-Born- 
Markovian master equation allows to determine the stationary long-time so- 
lution directly via straightforward diagonalization of the rate matrix. This 
allows to study the amplitude as well as the phase response for a wide range 
of parameters. Most importantly, we find pronounced resonances as well as an- 
tiresonances which are associated to multi-photon transitions in the resonator. 
We have found an interesting non-monotonous behavior of the antiresonance 
including a cross-over to a true resonant peak. This cross-over can be enhanced 
by the presence of the bath but is already inherent in the underlying coher- 
ent driven quantum system. Furthermore, we have found a clear separation of 
time scales in the dynamics how the globally stable state is approached start- 
ing from the metastable state. This tunneling process can be characterized 
by a single tunneling rate T which follows straightforwardly as the smallest 
non-zero eigenvalue of the rate matrix. If the control parameter being the 
frequency u is varied, resonant tunneling between the two states clearly is 
discerned as peaks in the tunneling rate. We hope that these rich features of 
the quantum Duffing oscillator will be found in future experiments. 
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Fig. 1. Amplitude A of the expectation value (x(t)) in units of 
x ° = \J mtn f° r var y m g the driving frequency u. Parameters are 
ksT = 0.lhujQ,a = O.lao,/ = 0.1/o,7 = 0.005u;o- Dashed line: Results of 
the classical Duffing oscillator at T = with the remaining parameters being the 
same. Inset: Amplitude A of the classical Duffing oscillator for varying driving 
frequencies. The driving strength / increases from bottom to top. 
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Fig. 2. Average energies E a (top), quasienergy levels e a: k (middle) and amplitude 
A of the fundamental mode for varying driving frequencies to. The remaining pa- 
rameters are a = O.lao, / = 0.1/o and 7 = 0.005a;o- 
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Fig. 3. (Color online) Response curve A(uo) for different values of the nonlinearity 
a. Remaining parameters are ksT = O.lfojQ, / = 0.1/o,7 = 0.005iuo. 
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Fig. 4. (Color online) Five-and six-photon resonance of the response curve A(uj) for 
different values of the driving amplitude /. Inset: Local extremum of the N— photon 
resonance for varying driving strength. The six-photon resonance develops from an 
antiresonance for small driving strengths to a true resonance for larger values of /. 
Remaining parameters are ksT = 0.lhujQ,a = 0.1ao,7 = 0.005il>o- 




17 



1 




Fig. 5. (Color online) Five-and six-photon resonance of the response curve A(u>) 
for different values of the temperature T. Inset: Local extremum of the N— photon 
resonance for varying temperature. The six-photon resonance develops from an an- 
tiresonance for low T to a true resonance for higher T. Remaining parameters are 
/ = 0.1/o, a = 0.1a , 7 = 0.005w . 




Fig. 6. (Color online) Five-and six-photon resonance of the response curve A(u) 
for different values of the temperature T. Inset: Local extremum of the N— photon 
resonance for varying temperature. The six-photon resonance develops from an an- 
tiresonance for low T to a true resonance for higher T. Remaining parameters are 
/ = O.l/o, a = 0.1a , 7 = 0.005w - 
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Fig. 7. (Color online) Phase shift <p(uj) for the five-and six-photon reso- 
nance for different driving amplitudes /. The dashed vertical lines mark the 
N = 5— and the N = 6— photon transition. Remaining parameters are 
T = 0.1T , a = O.Iqq, 7 = 0.005w - 
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Fig. 8. (Color online) Phase shift <p(u)) for different values of the temperature T. 
The dashed vertical lines mark the iV = 5— and the N = 6— photon transition. 
Remaining parameters are / = O.l/o, a = O.lao, 7 = 0.005wo- 
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Fig. 9. (Color online) Phase shift (p(uj) for different values of the damping constant 
7. The dashed vertical lines mark the N = 5— and the N = 6— photon transition. 
Remaining parameters are / = 0.1/q, a = O.lao, T = 0.1T . 
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Fig. 10. (Color online) Time-resolved dynamics of A for u = 1.167u;o (black solid 
line). Fast transient oscillations occur as "intrawell" relaxation in the metastable 
well. The long-time dynamics is governed by a slow exponential decay to the globally 
stable state characterized by a rate T (tunneling rate) . The red solid line shows a fit 
to an exponential e~ Tt . Here, a = O.lao, / = 0.1/o, fceT = O.I/kjq and 7 = 0.005u;o. 
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Fig. 11. (Color online) Tunneling rate T for the slow dynamics of the ampli- 
tude A for approaching the steady state for two different damping strengths. The 
peaks correspond to resonant tunneling in the dynamic bistability. Note that tun- 
neling is increased for increasing 7 (bath- assisted resonant tunneling). Moreover, 
a = O.lao, / = 0.1/o and ksT = OAHloq. 
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